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Abstract. We consider a two-dimensional convection model augmented with the rota- 
tional Coriolis forcing, Ut + U ■ WxU = 2kU'^, with a fixed 2k being the inverse Rossby 
number. We ask whether the action of dispersive rotational forcing alone, J7^, pre- 
vents the generic finite time breakdown of the free nonlinear convection. The answer 
provided in this work is a conditional yes. Namely, we show that the rotating Euler 
equations admit global smooth solutions for a subset of generic initial configurations. 
With other configurations, however, finite time breakdown of solutions may and actually 
does occur. Thus, global regularity depends on whether the initial configuration crosses 
an intrinsic, 0(1) critical threshold, which is quantified in terms of the initial vortic- 
ity, ojQ = W X Uq, and the initial spectral gap associated with the 2x2 initial velocity 
gradient, r/o := ^2(0) — Ai(0), Xj{0) = Xj{\/Uo). Specifically, global regularity of the 
rotational Euler equation is ensured if and only if 4:ku!o{a) + r]Q{a) < 4fc^, Va G IR^ . 

We also prove that the velocity field remains smooth if and only if it is periodic. An 
equivalent Lagrangian formulation reconfirms the critical threshold and shows a global 
periodicity of velocity field as well as the associated particle orbits. Moreover, we observe 
yet another remarkable periodic behavior exhibited by the gradient of the velocity field. 
The spectral dynamics of the Eulerian formulation, ^HIj reveals that the vorticity and the 
eigenvalues (and hence the divergence) of the flow evolve with their own path-dependent 
period. We conclude with a kinetic formulation of the rotating Euler equation. 

Key Words: Shallow-water equations, Rotational Coriolis forces, Critical Thresholds, 
Spectral gap, Kinetic formulation. 

AMS subject classification: Primary 35Q35; Secondary 35B30 



1. Introduction and statement of main results 

Finite-time breakdown is a familiar trademark of nonlinear convection mechanism. 
Consider the canonical example of an A^-dimensional system of free transport equations, 
Ut + U ■ VxU = 0. It follows — consult Corollarv 12.21 below, that the solution U{t,-) 
will lose its initial regularity at a finite-time if and only if an eigenvalue of the initial 
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velocity gradient crosses the negative real axis, i.e., iff there exists at least one eigenvalue, 
A(0, x) := A(^^^^p^), such that A(0,x) G IR~. Consequently, finite-time breakdown is 
a generic phenomenon for the free nonlinear transport. Thus, for example, irrotational 
initial data V^^ x f/(0, x) =0 — where all eigenvalues Aj(t, x) remain real, will necessarily 
lead to finite-time breakdown, except for non-generic cases where Xj{0,x) > 0, Wj,x, 
requiring, in particular, that the initial divergence is globally positive, ■ U (0, x) > 0. 
This general A^-dimensional scenario is completely analogous to the ID inviscid Burgers' 
equation, Ut + UUx = 0, where solutions of the latter will necessarily reach a finite-time 
breakdown except for the non-generic case of monotonically increasing initial data. 

Physically relevant models are governed by the fundamental Eulerian convection equa- 
tion augmented by proper forcing F, 

(1.1) Ut + U-VM = F. 

Here, there is a competition between the finite-time breakdown dynamics driven by non- 
linear convection and the balancing act of nonlinear forcing, F. Different models show up 
in different contexts dictated by the different modeling of such forcing. Three prototypes 
are dissipation, relaxation and dispersion. It is well known that if (jl.ll) is augmented 
with a sufficiently large amount of either dissipation or relaxation, then (jl.ip admits a 
global smooth solution for a rich enough class of initial data. In both cases of dissipa- 
tion and relaxation, global existence is secured by enforcing a sufficiently large amount 
of energy decay. Dispersive forcing, however, is different. The dispersive KdV equation, 
for example, Ut + UUx = Uxxx is a case in point. It admits global smooth solution while 
keeping the L^-energy invariant in time. In this paper, we study the regularity of the 2D 
convection model augmented by rotational forcing, 

(1.2) Ut + U ■ VxU = 2kJU, J : = 

subject to initial conditions, U{0,x) = Uo{x). Here 2k = e^^ where e is the Rossby 
number j^, e = expressed in terms of the characteristic length L, characteristic 
speed, U, and the amplitude of angular velocity Q of the rotating body, consult [HI . 
With these parameters the system evolves on a characteristic time scale t ~ L/f/. 

The system admits a global energy invariant in time, which is independent of the 
amplitude of rotation encoded by the constant k on the RHS of ()1.2j) . To see this, we 
note that ()1.2|) is formally equivalent to the extended 3x3 system, 

(1.3) dtp + Vx- (pU) =0, xGlR^ tGlR+, 

(1.4) dtipU) + Vx ■ {pU ®U) = 2kpJU, 

which are the usual statements of conservation of mass and Newton's second law, gov- 
erning the local density p = p(t,x) and the velocity field U := {u,v)(t,x). The usual 
manipulation, — l/2|[/px p.3|) + [/^ x p.4|) and the skew-symmetry form induced by the 
rotational forcing imply 



9t{lp\U\') + V. ■ ilpU\U\') = 2kp{U, JU) = 0. 
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The global invariance of the energy follows 

E{t) ■■=\ j p{t,x)\U{t,x)\^dx = E{0). 

The system (ll.3p - (jl.4j) coincides with a simplified version of the 2D shallow- water equa- 
tions, where the additional pressure terms are ignored. The only remaining forcing is 
the rotational Coriolis forcing, and our main quest in this paper is whether the action of 
dispersive rotational forcing alone prevents the generic finite time breakdown of nonlinear 
convection. The answer, outlined in Section 0] is a conditional yes. Namely, we show 
that ()1.2|) admits global smooth solutions for a subset of generic initial configurations, 
f/o. With other initial configurations, however, the finite time breakdown of solutions 
may - and actually does occur. Thus, global regularity depends on whether the initial 
configuration crosses an intrinsic, 0{\) critical threshold, which is quantified in terms of 
the initial vorticity, cuq := V^^ x [/q and the initial spectral gap, Fq := (A2(0) — Ai(0))^. 

Theorem 1.1 (Critical threshold for rotation forcing). Consider the 2D rotational flow 
m.^} with k > 0. Then the solution of / li.jjj) with initial data U{x, 0) = Uo{<y) remains 

x=a 

smooth for all time, — oo < t < oo, if and only if the initial data Uq satisfy 

(1.5) io{a) := 4:k[k-uJoia)]-To{a) > 0, Va G IR^ 

Moreover, as the smooth solution evolves along its particle path, then the vorticity, uj{t) = 
uj(t,a) and the eigenvalues Xj (and hence the divergence, d(t) = diVxU(t,a)) form a 
periodic orbit in phase space, with a path- dependent period, T = T{a), given by 

, , - 2 r^^ de 

(1.6) T =- 



kJ-n/2 {e,' + eo) + {eo'-eo)sine' 

Here 9q = 6{a) < 1 is determined by the initial data 
(1-7) Oq = —, po 



Several remarks are in order. 

1. We note that the critical threshold phenomena is independent of the initial divergence 
do := div^f/o. 

2. Let us point out that system (jl.2p could be viewed as a crossroad between the 2D 
shallow- water equations and the so-called 2D pressureless equations, e.g. [H 01131 EllSl ITB] . 
corresponding to Theorem 1 1 . II wit h A; = 0, 

(1.8) dtp + V-^- (pU) =0, xGlR^ tGlR+, 

(1.9) dt{pU) + V,-{pU®U) = 0. 

According to Theorem II. H the pressureless system admits a global smooth solution for- 
ward (respectively, reversible) in time, if and only if Xj{0) ^ IR~, respectively iff Xj{0) ^ IR. 
The latter is equivalently expressed in Theorem II. II requiring Tq := (A2(0) — Ai(0))^ < 0. 

3. In particular, (jl.2j) does admit global smooth solutions with negative initial diver- 
gence in contrast to the free transport {k = 0) equation discussed in the introduction. 
It follows that rotation prevents finite time breakdown, either by a large Coriolis forcing 
(/c >> 1) or a large initial rotation (Fq << 0). 
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4. If we set y-independent initial data, then ()1.2j) is reduced to the one-dimensional 
system 

Ut + uux = 2kv, 
Vt + uVx = —2ku, 

with critical threshold ('UQ(a))^ — Akv^la) < Ak"^. To interpret Theorem 11.11 in this sim- 
plified setting, we observe that the gradient {uj,d) := {—Vx,Ux) solves a coupled system 

(1.10) {dt + udx)uj + du = 2kd, 

(1.11) {dt + udx)d + d^ = -2kuj, 

and a straightforward computation reveals the path-dependent invariant along the particle 
path, X{t) = u{X,t),X{0) = a. 



d"^ + uj^ ' ' dl{a) + LJo{ay 

The critical threshold statement in this case reads Bo > 1, stating that the gradient 
{uj, d) forms a closed elliptical orbit in the phase plane (whereas Bq < 1 corresponds 
to unbounded parabolic/hyperbolic orbits). Following the analysis in Sectional we also 
obtain a path-dependent period for the gradient 



k (%' + ^o) + (%' - Oo)stne' Vm+l 

Such path-dependent period of the gradient refiects the fact that its governing system 
fll.lO|) . ()l.lH) . is a nonlinear perturbation of the harmonic oscillator. As the Rossby number 
approaches zero, however. A; >> 1, 6*0 ~ 1, and the above path-dependent period T is 
approaching the global inertial period n/k (the harmonic oscillator period). 

5. As we shall see in Theorem II . 21 below, sub-critical initial data yield smooth velocity 
fields, U{X{t)), with time period T = it/ k, or — expressed in terms of the original non- 
scaled time units, a period T = T U / L = rr/Q. This period of the particle orbits is related 
to the global revolution of the plane. Theorem 11.11 points out yet another remarkable 
property for a portion of the gradient of U(t, ■), namely, the divergence d(t,a), and the 
vorticity uj(t, a) which exhibit a local, path-dependent period dictated by the unique 
initial parameter, Sfcpo- It is instructive to compute the period predicted in Theorem ll.il 
using configurations similar to those encountered in various applications. Let us illustrate 
a couple of examples taken from [IT. For the Gulf Stream, with the Rossby number 
e = 0.07, L = 100 km and U = Im/sec, we find that the vorticity and divergence of the 
fiow keep repeating themselves every T = T L/U ~ ll.Yhrs; for the weather system we 
have e = 0.14 with L = 1000km, U = 20m/ sec, and the vorticity/divergence exhibit a 
period of T = T L/U ~ 12.2hrs. It is also interesting to see how this path-dependent 
gradient period be influenced by the small Rossby numbers. After rescaling we may 
assume initial configuration such that do ~ ^^o ~ 1? for which a small Rossby number 
yields tQ ~ Ak"^, Pq ~ 2k and hence 6q ~ 1. Restored in terms of the original time scale, 
T = T L/U, the period is given by the 6'o-dependent elliptic integral, 

2_ r^' de TT 



n y_,/2 (^-1 + ^o) + (^0^^ - eo)stne 
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which is close to the inertial period when the Rossby number is small. For the earth core, 
for example, we have T = 11.95/zrs with e = 2 x 10~^, (L = 3000/cm and U = 0.1cm/ sec), 
whereas for Jupiter's Red Spot we have a Rossby number e = 0.015 (with L = lO^km 
and U = 0.1cm/ sec) and the velocity gradient period T ~ 5.13hrs. We should point out 
the difference between the period of the velocity field vs. the velocity gradient periods, 
which is due to our tracking of the flow dynamics along the particle path. 

Of course, one should not expect the current cartoon model to provide a faithful descrip- 
tion of the full model, as other forces ignored at this stage — magnetic forces, pressure, 
etc, will play the important role. Nevertheless, the above rough approximations are in- 
teresting for their own sake, in particular, since the flow is predicted to be periodic once 
smooth solutions are secured for subcritical initial data. Surprisingly, the periods com- 
puted above fall within the physical range. It will be challenging to refine the estimate, 
by taking into account the other forces which augment with the rotation model. 

To put our study in a proper perspective we recall that there has been a considerable 
amount of literature available on the global behavior of nonlinear convection driven by 
rotational forcing and related problems, from rotating shallow-water model ^21 ^1 122] 
to rotating incompressible Euler and Navier-Stokes equations j3 El El EI- The common 
feature of the flows studied in this context are rotation dominated flow for which the 
Rossby number e is small. It is well known that large-scale atmospheric (or oceanic) fields 
are in permanent process of Rossby (or geostrophic) adjustment [20]. The flow structure 
has been extensively studied in terms of e, say in 112], based on the averaging of the 
interaction of the fast waves of the rotating Euler equation, two dimensional structures 
were shown to emerge in the limit e ^ 0; for bounds for the vertical gradients of the 
Lagrangian displacement that vanish linearly with the maximal local Rossby number jU]; 
for a nonlinear theory of geostrophic adjustment for the rotating shallow-water model 
under the assumption of the smallness of the Rossby number [22]; consult [12] for the 
analysis of an approximation of the rotating shallow-water equation. 

When dealing with the questions of time regularity for Eulerian dynamics without 
damping, one encounters several limitations with the classical stability analysis. Among 
other issues, we mention that 

(i) the stability analysis does not tell us how large perturbations are allowed before 
losing stability - indeed, the smallness of the initial perturbation is essential to make the 
energy method work, e.g. the 3D Navier-Stokes equation (IGj : 




(ii) the steady solution may be only conditionally stable due to the weak dissipation in 
the system, say in the ID Euler- Poisson equations [TT| . 

To address these difficulties we advocated, in |1H[T71 ITH]. a new notion of critical thresh- 
old (CT) which describes conditional stability, where the answer to the question of global 
vs local existence depends on whether the initial configuration crosses an intrinsic, 0{1) 
critical threshold. Little or no attention has been paid to this remarkable phenomena, and 
our goal is to bridge the gap of previous studies on the behavior of rotational Euler equa- 
tions, a gap between the regularity of Eulerian solutions in the small and their finite-time 
breakdown in the large. The critical threshold (CT) was completely characterized for the 
ID Euler-Poisson system in terms of the relative size of the initial velocity slope and the 
initial density; consult f?! |21] for the CT for the convolution model for conservation laws; 
Moving to the multi-D setup, one has first to identify the proper quantities which govern 
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the critical threshold phenomena. In JHl we have shown that these quantities depend in 
an essential manner on the eigenvalues of the velocity gradient matrix, XCVxU). 

The critical threshold for the current rotation model can be also obtained, in a straight- 
forward manner, through a Lagrangian flow formulation. This is summarized in the Theo- 
rem 1.2 below. We should point out that it was the spectral dynamics analysis of XCV^U) 
that led us to the CT formulation in the first place, which in turn was then sought within 
Lagrangian formulation. In Section El we prove 

Theorem 1.2 (Flow map for rotation forcing). The flow map, Xa '■= = U{Xa), 
-^q(O) = a associated with M.^) . = 2k J X^ is given by 

Xa{t) = j^J-'e^''''Uo{a) + a- ^J-'Uo{a). 

For sub-critical initial data, m.5]} . this flow map is invertible and periodic with an inertial 
period T = ir/k. The velocity field U{Xc,{t)) = Uo{a) + 2kJ{Xa{t) — a) shares the same 
inertial period. 

At this point, one may wonder whether this inertial period is none other than the planet 
rotation. Actually the two are not the same; the rotating plane completes one revolution 
in a time equal to while the period of the particle path expressed in the original non- 
scaled variables is T = TL/U = ^. Thus, the particle goes around its orbit twice as the 
plane accomplishes a single revolution, which is consistent with the observation in pH] . 

Finally we conclude in Section IHl with a kinetic formulation of the current rotation 
model. 

Theorem 1.3. The rotation model il.^) admits for the following kinetic formulation 

dtf + i- + 2kJi . V^/ = ^(M - /), 

where My^u}{C) is the Maxwellian given by 

where p and U are macroscopic density and velocity respectively, and T being an arbitrary 
fixed temperature. 

In Section|3]and|Slbelow, we quantify the same critical threshold using the Eulerian and 
Lagrangian formulations, and it would be of interest to derive the same critical threshold 
directly using the kinetic formulation in Theorem 11.31 

2. Spectral dynamics 

We consider a general nonlinear transport equation (jl.lj) . Ut + U ■ V^U = F, and we 
trace the evolution of VxU in terms of its eigenvalues, A := X{'VxU)(t,x). The following 
result is in the heart of matter. 

Lemma 2.1. [ Spectral dynamics, Lemma 3.1] ] . Let A := A(V^f/)(t, x) denote 
an eigenvalue of VxU with corresponding left and right normalized eigenpair, {i,r) = 1. 
Then A is governed by the forced Riccati equation 

dtX + U- VxX + X^ = {£, VxFr). 
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As an immediate corollary we obtain the precise description for finite time breakdown 
of free nonlinear transport. 

Corollary 2.2. [Finite time breakdown of free transport, |^ Lemma 4.1]]. The free 
nonlinear N- dimensional transport 

(2.1) dtU + U-V^U = 0, xelR^, 

admits global smooth solution forward in time, t > 0, if and only if the eigenvalues of its 
initial velocity gradient, X := \(VxU), satisfy \{0,x) ^ IR^. Likewise, it admits a globally 
smooth, time-reversible solution for — oo < t < oo iff X{0,x) ^ IR. 

For the proof, we note that the eigenvalues, governed by the homogeneous Riccati 
equations, propagate along the particle path x = x(t, a), 

X{t,x) = -^^^. 

We note in passing that the rotational system ()1.2|) is to the full shallow-water equations 
as the free transport model ()2.1|) is to the full Euler equations. The existence of a critical 
threshold phenomena associated with global linear forcing model was first identified by us 
PH] ! although the exact configuration cannot be obtained in such generality. The current 
paper provides a precise description of the critical threshold for the 2D rotational system 
(II. 2|) . In particular, we use the Spectral Dynamics Lemma to obtain remarkable explicit 
formulae for the critical threshold surface summarized in the main Theorems 11.11 Taking 
the gradient of the velocity equation ()1.2|1 . we find that the velocity gradient field, VxU, 
solves the following matrix equation 

(2.2) dtiVxU) + U ■ V,.(V.t/) + (VxUf = 2kJ\/xU. 

Using the spectral dynamics Lemma f2. 11 we obtain the spectral dynamics equations 

(2.3) dtXi + U ■ Vx.Ai + Xl = 2kXi{k, Jr.), 

(2.4) dtX2 + U ■ VxX2 + Xl = 2kX2{l2, Jrs), 

where Xi,i = 1,2 are eigenvalues of the velocity gradient field VxU associated with left 
(row) eigenvectors k and right (column) eigenvectors rj. Since J is skew-symmetric we 
have Jri = ailj and Jr2 = a2lj ■ Noting that l2r2 = hri = 1, one then has ai = (r2, Jri) 
and 

"2 = {ri, Jrs) = -(ra, Jri) = -ai. 
Therefore, (|2.3p - (j2.4j) now read 

(2.5) dtXi + U-VxXi + Xl = 2kXi{r2,.M){h,l2), 

(2.6) dtX2 + U ■ \/xX2 + Xl = -2kX2{r2, Jr,){h, k), 

from which we deduce that the spectral gap rj := X2 — Ai and divergence d := X2 + Ai, 
satisfy 

dtr] + U- VxV + dr] = -2kd{r2, Jri){h, h) 

and 

(2.7) dtd + U- Vxd + = -2kv{r2, Jri) {h, h). 
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On the other hand, differentiation of ()1.2|) yields the V^^fZ-equation ()2.2p . i.e., 

(2.8) id,+u.v.) (i^ iy) + ( "^+;-"^ ^ j^^i^, ) =2^( -i 

\ '-'y J \ "'f^a; '^x'^y ' / \ "'x '^y 

which in turn - using the LHS of ()2.5|) - ()2.fij) to express \l + Xl = {cP + rf)/2, leads to 

(2.9) dtd+U =-2fccj, 

(2. 10) 9fCj + f/ ■ V:,.cj + c/cj = 2A;d. 
Equating the expressions on the right of ()2.9|) and ()2.7|) we find 

-2ku = -2kr]{r2, Jri){lij2)- 

Thus, the scaled product of the eigenvectors measures the ratio of vorticity over the 
spectral gap in following manner 

(2.11) {r2,Jn){h,l2) = -. 

V 

When the spectral gap r] shrinks to zero, the scaled product becomes unbounded due to 
the degeneracy of eigenvectors. When the vorticity u shrinks to zero, (j2.1H) recovers the 
symmetry of VxU which is reflected through the orthogonality of ii and £2- Equipped 
with the above relations we come up with a closed system for [uj, d, rj) along the particle 
path (here and below ' = dt + U ■ V^) 

u' + du = 2kd, 
2 

r] + di] = —2k — . 

V 

Note that the spectral gap may become purely imaginary when eigenvalues are complex. 
To avoid the discussion on the complex solution of the above system, we introduce the 
following real variable 

r := r]^ 

Using the above equations we have 

r = 2r]r]' = 2d[-2ku - T]. 

Note that the sign of 2k indicates the direction of the rotational forcing, and the vorticity 
measures the rotation in the flow. In order to combine these two effects we introduce 
{p := Ak"^ — 2kuj, and thus obtain a closed system for W := (y?, d, T)'^ 

(2.12) Lf' = -dip, 

(2.13) rf' = _^JLL + ^_4fc2, 

2 



(2.14) r = 2d[^-4F-r]. 
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We shall use this system to describe the dynamics of the velocity gradient field. Lin- 
earization of the above system around W* = {ip* ,d* ,T*)~^ gives the linear system W = 
A{W*){W -W*) with 

/ -d* -if* \ 
A = I 1 -d* -i 

\ 2d* 2ip* -2d* J 

The corresponding eigenvalues of A at critical points 0,r*) are Ai = 0, A2,3 = 
±y/—2ip*. The classical stability analysis based on linearization is not sufficient to predict 
the global time dynamics. 



3. Material Invariants 

It follows from equations ()2.12|) and ()2.14|) we obtain 

dip —if) 

dT ^ 2((^-4A;2 -r)' 

which upon integration gives the first material invariant 

(3.1) =Co, Co = Co{a):= ^ . 

This material invariant enables us to reduce the full system ()2.12|) - ()2.14p to the following 

system 

(3.2) if' = -ifid, 

(3.3) d' = -^[d^ + - Coip% 

In order to have global bounded solution it is necessary to assume Co(a) > 0, i.e., 

(3.4) To < 2(^o-4A;^ =4fc(fc-cJo), 

for otherwise, ()3.3j) will be majored by the Riccati equation d' < — ^ — 2A;^, which would 
lead to the finite time breakdown. As we shall see below in §4, the positivity of Co{a) is 
also sufficient for the global bounded solution. 

Another material invariant is obtained in the following manner. Following Lemma 
2.2], we set g = d"^ which yields 

dq d' q + Ak^-Coip' 
—— = 2d — = . 

dip if' if 

Integration yields a second material invariant 

and together with ()3.1|) we end up with a second independent material invariant 

d^-T + 2ip 



In summary we have 



Vo{a). 
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Lemma 3.1. Let ip := Ak'^—2kuJ, d := divU = tr{WxU) andV := (A2 — Ai)^ he the solution 
of the dynamical system i2.1^) - ^F^ ), associated with the rotational system il.^) . Then 



we have the following material invariants along particle path {t,Xa{t)), 
(3.6) 



(3.7) 



d'^-T 



_ dl{a) - To{a) 
(t,x^{t)) (/)o(a) 



4. Critical Thresholds 



As we observed earlier, the positivity of condition, Co(a) > 0, is necessary for global 
bounded solution, for otherwise 

d> < -l[Ak'' + d% 

which would imply that c?, and hence yj, become unbounded in a finite time. We shall show 
that the same positivity condition, Co (a) > 0, is in fact also sufficient for the existence 
of global bounded solution. For Co > 0, the reduced system ()3.2|l - ()3.3j] has two unique 
equilibrium points in the finite plane, [ip*^,d) = (±^^,0). The local behavior of the 
solution depends on the properties of these critical points. We note that since ip = 
is an invariant set, then (pQ(p{t) > for all time, and we therefore concentrate on the 
solution behavior for yjo > 0, with the other case of <y9o < being handled similarly. On 
the right plane > 0, the coefficient matrix of linearized system of ()3.2|l - ()3.3j) around the 
equilibrium point (v5+,(i) = (-^,0) is 



-y,* 

with purely imaginary eigenvalues, ±(^/CoV^+)^- This means that the bounded trajectory 
is possibly a periodic solution or limit circle. Observe that if {(p{t),d{t)) is a solution to 
()3.2|) . ()3.3|) . so is {{p{—t),~d{—t)). Such symmetry implies that (</?+, 0) is a center and 
the trajectory in a neighborhood of this equilibrium point is periodic. 
In order to clarify the global behavior of the flow around the center, we appeal to the 
material invariant ()3.5|1 which we rewrite as 

d' + iVCoip - 2ky 

V+{(p,d) := . 

V^{-) is a positive deflnite Liapunov function for > 0, and achieves its global minimum, 
= 0, at the equihbrium point (93I = ^=,0). A family of closed orbits in the phase 



plane {(p, d) can be expressed as the level set curve 

V+{({), d) = Const > 0, 

since V+ is material invariant in the sense that = 0. Similarly, on the left plane < 0, 
one may use the Liapunov functional 

d^ + iVCo^ + 2kf 
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whose level set curve determines a family of closed orbit on the left plane centered around 
{lp*_ = 0). The global behavior of the solutions is summarized in 

Lemma 4.1. [Bounded solutions are periodic] 

{i} Given initial data {ipo, do, Tq) for the system \2.12^) - \2.1Ji\l , it admits a global hounded 
solution if and only if the initial data lie in the suh- critical region where 

(4.1) To < 2(^0 - 4A;l 

{ii\ Bounded solutions of I12.1^} - ^2.14\ ) ^'^c necessarily periodic. The periodic orbit on the 
right plane (f > lies on the ellipse d"^ + (-s/CqV' ~ 2A;)^ = {Vq — Ak^/CQ)ip, with Co and Vq 
are determined by the initial data, 4,'/. i)) and ( 1,7. 

Proof. By ()4.1|) . the initial data for the bounded orbit must satisfy Fq < 2ipQ — 4/c^, i.e. 
Co > 0. From the fact that V{ip, d) = Const > is the material invariant, we find that F 
can be expressed as 

d' + {^oV - 2kY = {Do - 4k^o)y^- 
This is an ellipse, along which the corresponding solution (</? = 4fc^ — Akuj, d) being 
periodic. Moreover, the invariants ()3.6p . ()3.7|) imply that F shares the same period along 
the path. It follows that the boundedness of the divergence implies the boundedness of 
the whole velocity gradient |19j [Lemma 2.1]: the anti-trace of ()2.8p r := + Uy satisfies 
r' + rd = —2ks followed hj a. s := — f^-equation, s' + sd = 2kr yield the boundedness 
of W. In fact solving the above two equations in terms of the divergence d we obtain 

- = tan(tan~^(ro/so) - 2kt), r'^ + s^ = (rg + si) exp ^-2 d{^)d^^ . 

By the periodicity of d and its symmetry about the axis d = Owe see that s^+r^ shares the 
same period with d; This is also clear from the identity + = F + tu^. A combination 
of the above facts lead to 



r = sm (tan ^(ro/so) — 2kt)VT + lu"^, s = cos(tan ^(tq/so) — 2kt)VT + uS^. 

As a product of two periodic functions with period vr/fc and T respectively, the anti-trace 
r, also the anti-vorticity s, is periodic only if the ratio of these two periods, i.e., 

"1 di 



(0o + ^^o"') + (^o-^o"')cos(7rO 

is a rational number. In this case the overall gradient V^^f/ is also periodic with period 
niTt/k for some integer m. This completes the proof. □ 

Once we identified bounded solutions as periodic, the next step is to seek the period 
for each periodic orbit. 

Lemma 4.2. The period of each bounded orbit associated with \2.1^) - ^.14}) is given by 

7^ 2 r/^ dO 

(4.2) T = - 



9o + 9^^ + {e^^-9o)sine' 
Here Oq = 9o{a) < 1 is given by 



(4.3) 9„ := 

VI + Skp„ + 1 
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where Pq is determined by the initial data 



Po[a) 



Proof. Due to the symmetry it suffices to compute the half period. The intersection points 
of the elhpse V{(f, d) = V{(f±, 0) with d = can be written exphcitly in terms of the initial 
data 

(4.4) ^- = ^^0, V^ = ^0o'- 

Using the equation 

dip 

Tt = 

where along the trajectory from 0) to 0) in the lower-half phase plane [ip, d) we 
have 

d = -VDqlp - Cov?2 _ 4k'^ = -a/C^a/ -ip){ip- if^). 

Therefore 

d^ 2 /"^+ ds 



(4.5) T = 2 , - . 

-ipd VCq Jy,_ s^/{lp+ - s){s - ip-) 

Let s = "^-^^^ + f+^E^T^ and using the expression oi ip- and in ()4.4|) one has 

- 4 /"^ dr 
T 



/Co 7-1 [v5- + + (V5+ - ip-)T]^/l - 
4 /""/^ 



+ + (^+ - 'P-)sine 

2 r'"^ de 



kJ-./2 eo + e^' + {e,'-eo)szne' 

This gives the desired result. □ 

5. Flow Map 

For the smooth flow we may further study the structure of the flow map. Assume 
X = Xa{t) is the flow map started at the initial position a, then one has 

dX 

X^:=^ = UiX^), X„(0) = a, 
and momentum equation can be written as 

= 2kJXa- 

Integration once gives 

Xa = Uo{a) + 2kJ{Xa - a), 

where Uo{a) is the initial velocity at location a. The above equation leads to the flow 
map expression 

Xa{t) = ^J-\^^''U^{a) + « - ^ J-^f/o(«). 
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This flow map determines the unique smooth velocity field if and only if the indicator 
matrix, 

remains nonsingular. Noting that J ^ = — J and 

cos{2kt) sin{2kt) 
—sin{2kt) cos{2kt) 



we have 



J f 1 - cos{2kt) -sin{2kt) , 

2k \ sm{2kt) 1 - cos{2kt) ' 



Hence, with U = {u, v) we find 

^t-r -9hT ^( ^^^(2^^) 1 - cos{2kt) \^ rj 
2kV -2kl+y - 1 sin{2kt) ) 

2k + vqx + UxSin{2kt) — voxCos{2kt) Voy + uoySin{2kt) — voyCos{2kt) 
—uqx + voxSin{2kt) + Uoxcos{2kt) k — uoy + VQySin{2kt) + uoycos{2kt) 

A careful calculation gives its determinant as 

det{2kT) = 4k^ - 2kuo + 2det{V JJq) + {2kuJo - 2det{V JjQ))cos{2kt) + {2kd)sin{2kt) 

Thus r(t) remains nonsingular for all time if and only if det{2kT) ^ 0, i.e., 

4fc2 - 2kuQ + 2det{VM 



^ \^-^{2kuJo - 2det{VUo)y + 4k^dl {2kojQ - 2det{V + Ak^dl 
which is equivalent to 

(5.1) (4A;2 - 2kujQ + 2det{V > C^kcuo - 2det{V JjQ)f + Ak'^dl. 

The spectral gap Fq relates the determinant and the divergence as 

To = dl-Adet{VJJ^). 
This when applied to the above inequality (j^.lj) gives 

To < Ak^ - AkiUo, 

which is exactly the critical threshold in phase space stated in Theorem ll.il 

6. Kinetic formulation 

This section describes a kinetic formulation for the rotational model (jl.2j) in terms of 
a density function, / = f{t,x,^) governed by the BGK model 

(6.1) dtf + ^-Vxf + 2ke-V^f = -^{M-f), ^^:=Je- 

Here M = M^p^uy{^) is the Maxwellian given by 

M = ^e-l«-^l^/^, e=(ei,6)GlR'- 
V ttT 
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The fixed temperature T, plays no role in this pressureless model. The connection between 
the distribution function / and macroscopic flow variable is realized in terms of the usual 
moments of density p, momentum m = pU and total energy E = p|f/p/2, 

|2xT 



(p, pU, E) 



The conservation principle for mass, momentum and energy during the course of particle 
collisions requires the equilibrium to satisfy the compatibility condition 



while the rotational forcing is introduced through the potential 

j 2ke ■ V^/^^(Ocie = (0, -2kpJU, 0). 

Indeed, a straightforward computation yields / 2kS,~^ ■ V^fdC, = 0; for the momentum 
equation we compute 

J 2ke ■ V^fidi = 2k j ^(6%/ - ^AJ)di 



2k 



-pv 
pu 

2kpJU, 



while the presence of this forcing does not change the energy equation since 

2k j e^.Ve/|eprfe = 0. 
The first three moments of then yield the equivalent extended system of (|1.3p . (jl.4|) . 

dt 







pU 1 


+ Vx - \ Fm 


V Fe 




The corresponding macroscopic fluxes are 

{Ep,F^,FeV := j ^^(OeM, 
and under the closure / = M^p^uj we conclude 

(Fp, Fm, FeV = [pU, pU ® U, p^ 
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